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AN ANALYSIS OF THE RELAXATION OF LAMINAR BOUNDARY LAYER 
ON A FLAT PLATE AFTER PASSAGE OF AN INTERFACE WITH 
APPLICATION TO EXPANSION -TUBE FLOWS 

By Roop N. Gupta* 

Langley Research Center 

SUMMARY 

The relaxation of the accelerating-gas boundary layer to the test -gas boundary layer 
over a flat plate in an expansion tube is analyzed. Several combinations of test gas and 
acceleration gas are considered. The problem is treated in two conically similar limits: 
(1) when the time lag between the arrival of the shock and the interface at the leading edge 
of the plate is very large, and (2) when this lag is negligible. The time -dependent 
laminar-boundary-layer equations of a binary mixture of perfect gases are taken as the 
flow -governing equations. This coupled set of differential equations, written in terms of 
the Lam-Crocco variables, has been solved by a line -relaxation finite -difference tech- 
nique. The results presented include the Stanton number and the local skin-friction coef- 
ficient as functions of shock Mach number and the nondimens ional distance -time variable. 
The results indicate that more than 95 percent of the test -gas boundary layer exists over 
a length, measured from the leading edge of the plate, equal to about three-tenths of the 
distance traversed by the interface in the free stream. 

INTRODUCTION 

An expansion tube (ref. 1), a facility for producing high -enthalpy gas flows, is one 
of the aerodynamic testing devices in which the test -flow duration is very brief. Conse- 
quently, the successful use of such a device makes it necessary to understand the nature 
of the flow development over the test model and the time required to attain steady flow 
conditions, and, in turn, requires a detailed knowledge of the time -dependent viscous- 
flow processes. In the case of an expansion tube, the model is first immersed in the 
flow of the accelerating gas prior to the arrival of the test gas and, thus, it is necessary 
to know the time required for the accelerating -gas boundary layer to relax to the test- 
gas boundary layer. The inviscid flow relaxes more rapidly than the boundary layer. 

For the case of a flat plate, the inviscid flow may be assumed to be fully developed in the 
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(d) After secondary diaphragm burst, t = t-5. 



(e) During testing, t = t c . 

Figure 1.- Schematic diagram of an expansion-tube operation 
cycle. (Taken from ref. 1, circled numbers are explained 
therein. ) 






time taken by the particle in the flow to travel from the leading edge to the last down- 
stream point of interest on the plate. The boundary layer requires longer to become fully 
developed because of the viscous processes. Under certain conditions when the testing 
time is very short, the steady -state conditions may not be reached. In these cases, an 
evaluation of experimental data requires the prediction of shear stress, heat transfer, 
and other parameters in the unsteady flow region. 

In the present study, the shock-induced flow over a semi-infinite flat plate in an 
expansion tube is analyzed. A schematic representation of the operating cycle of an 
expansion tube is presented in figure 1. The apparatus is divided by two diaphragms into 
three sections. The driver section contains a gas at high pressure. The driven section 
contains the test gas, and the expansion section contains an accelerating gas at low pres- 
sure. The test model is located near the downstream end of the expansion section. The 
test region is the area labeled (IT) in figures 1(a) and 1(e). Reference 1 should be con- 
sulted for more details on the operation of an expansion tube. 

For quantitative evaluation of the boundary -layer quantities over the plate surface, 
nitrogen is considered to be in the test section, and helium and argon are considered to 
be in the accelerating section. The reason for using argon in the accelerating section is 
to obtain the effect of an accelerating gas with a larger molecular weight on the relaxa- 
tion of the boundary layer. 

The relaxation of the helium boundary layer and the argon boundary layer to the 
test -gas (nitrogen) boundary layer is analyzed by solving the time -dependent boundary- 
layer equations of species concentration, momentum, energy, and continuity, along with 
the equation of state for a binary gas mixture. The appropriate transport and thermo- 
dynamic properties for the perfect gas mixture are evaluated by employing the Monchick 
exponential repulsive potential. (See ref. 2.) The governing equations are treated in 
the transformed plane by using a conical similarity variable to obtain the self -similar 
solutions. The resulting set of coupled second-order nonlinear partial differential equa- 
tions are solved by a line -relaxation technique which is a method that has been used 
extensively for elliptic problems. The governing boundary -layer equations employed 
herein, although parabolic in nature, require boundary conditions at the leading edge as 
well as at the interface because of a change in sign of the coefficient of the convection 
term in this region. This type of behavior has been termed "singular parabolic" by 
Gevrey (ref. 3) and Lam (ref. 4). 


SYMBOLS 

A ratio of shock velocity to free -stream velocity, U s /U 0 

a fraction of model length 


3 



P MP 

' (MP)w 


°12’ 12 

molecular thermal diffusion factor parameters 

c3C i 


c i 

mass fraction of heavy gas , p^Jp 

c 2 

mass fraction of light gas, p^fP 

c i 

local skin -friction coefficient, 

2 v w U ° )n 2 

C P 

specific heat of mixture at constant pressure 

c v 

specific heat of mixture at constant volume 

»i T 

thermal diffusion coefficient of heavy gas 

D 12 

binary coefficient of diffusion 

F 

function 

G 

Blasius shear function, (p{0 ,0) 

H 

k - h w 

dimensionless enthalpy difference, , 

a w 

Hq 

, h o " h w 

dimensionless enthalpy difference behind the shock in free stream, ^ 

h g 

Blasius enthalpy function, H(0,/3) 

h m 

Mirels enthalpy function, H(1,0) 

h . 

local enthalpy of mixture; also step size in numerical method 

hi = c^i 

where i = 1, 2 

ho 

local enthalpy behind shock in free stream 
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J 

k 

M 

M 

M s 

m 

N L e 

Npr 

N St 

n 

P 

q 

R 

R’ 

**w,x 

T 

T 0 

T r 

T s t 


thermal diffusion factor, ^DjT also iteration parameter 


thermal conductivity of mixture 

Mirels shear function, 0(1, /3) 

molecular weight 

shock Mach number 

count of steps in a -direction 
pE> 12 c 

Lewis number, - 


c /J. 

Prandtl number, -2— 
’ k 


Stanton number, -q 


w/ 


"( P w)N2 Uo ( C p)N 2 f Tr ' Tw ^ 


count of steps in /3-direction 


pressure; also iteration parameter 


heat -transfer rate 


gas constant 

universal gas constant 

P w Uq x 


Reynolds number, 


M 


w 


temperature, K 

free-stream temperature in region between shock and interface 
recovery temperature 

stagnation temperature in free stream (having temperature TJ 



Tj temperature of gas in front of shock 

free-stream temperature in region between interface and expansion fan 

t time, measured so that at t = 0 the shock wave is located at leading 

edge of plate 

t’ time, measured from the moment secondary diaphragm is ruptured 

t* time, measured so that at t = 0 the interface is located at leading 

edge of plate 

t transformed time variable in Crocco plane, identical to t 

U Q free-stream velocity 

U s shock velocity 

u velocity component parallel to plate surface in boundary layer 


u s u(x,y,t) 

v velocity component normal to plate surface in boundary layer 

, ^ 9 In c 1 Dj 7 a in T 

V 1 ~ ' d 12 ay ' pcj 9y 

x distance from leading edge along plate in x,t coordinate system 

x . distance from secondary diaphragm along center line of expansion tube 

x * distance from leading edge along plate in x* ,t coordinate system 

x transformed distance variable in Crocco plane, identical to x or x* 


normal distance from plate 


a. 


a 


* 


conical coordinate in MR limit, 
conical coordinate in BL limit, 
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a ., = — 
crit K * 

o 

|3 dimensionless velocity, — 

u o 



y 

y 

_* 

y 
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P 

T 

0 

0 


^2 y l ' 1 M 1 

y 2 - 1 Vi m 2 


ratio of specific heats, 


dimensionless time in MR limit, defined by equation (A16) 


P w u o 2 t* 

dimensionless time in BL limit, 

Mw 


thickness of velocity boundary layer 


“ ln~ dy where e denotes edge of boundary layer 
u o/ p e 


dimensionless y-coordinate 


momentum thickness , \ (l - ~ Y-^ ~ dy 

J 0 V U o/ u oP e 


+ (1 - c) 


dynamic viscosity of mixture 
kinematic viscosity of mixture 
density of mixture 
shear stress 

shear function defined by equation (A25) 
dimensionless shear stress defined by equation (A21) 
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Subscripts: 


1 heavy gas 

2 light gas 

e edge of boundary layer 

max maximum 

N 2 evaluated for nitrogen 

w evaluated at wall 

Abbreviations: 

Ar argon 

BL Blasius 

He helium 

MR Mirels 

N 2 nitrogen 


BACKGROUND 

Several authors (refs. 4 to 19) have evaluated the various boundary -layer parameters 
for a shock tube having similar unsteady characteristics. References 5 to 14 are related 
to the problem of a boundary layer on the walls of a shock tube, whereas references 4 
and 15 to 19 discuss the boundary layer on a semi-infinite flat plate which develops behind 
an advancing shock wave. This latter problem is relevant to the present study and will be 
discussed in more detail. 

Trimpi ancTCohen (ref . ~6) used~an~integral-method-fer-calculating-the-laminar bound - 
ary layer for the entire flow in a shock tube. They integrated the boundary-layer equa- 
tions in the direction normal to the wall bounding the flow. These equations were further 
transformed into a conically self -similar coordinate system. The resulting hyperbolic 
differential equations were solved by the method of characteristics utilizing an integral 
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technique at the discontinuities (contact surface, and so forth) in order that the charac- 
teristic solution may proceed across the discontinuities. 

Mirels (ref. 8) and Ackroyd (ref. 14) used a shock-fixed coordinate system in which 
the boundary -layer flow is steady and similar solutions are valid. Although the shock- 
based coordinate system appears to be very appealing for the flow conditions which can 
be approximated mathematically by infinitely long walls, it does not facilitate the solution 
of the governing equations for the case of a shock moving over a sharp leading -edge plate. 
The solution obtained in the shock-fixed coordinate system for such a problem is valid 
only in the vicinity of the shock and cannot be expected to satisfy the boundary condition 
at the leading edge of the plate. 

Cohen (ref. 9) analyzed the boundary layer developed by a centered expansion fan 

advancing into a stationary fluid by reducing the boundary -layer equations to similarity 

form by using the conical similarity variable, 4 = 1- zr-r- He obtained the three -term 

U 0 t 

power-series expansions of the stream function and the enthalpy in terms of 4 (the 
distance behind the wave head) which are valid for small values of 4 only. Becker 
(ref. 10), using a numerical continuation procedure, was successful in extending Cohen's 
results for velocity and temperature profiles up to 4 = 0-9. Basically, Becker’s continua- 
tion procedure requires taking an initial profile given at some 4* with which the new 
profile is computed at 4 > 4 + . However, at the tail of the expansion wave 4 0 where the 
derivatives of the free -stream quantities are discontinuous, the continued profile fails to 
satisfy the wall compatibility conditions for 4 > I n order to overcome this difficulty, 
Becker added correction terms having proper asymptotic behavior far from the wall and 
approximating the correct result near the wave tail 4 0 * Uniqueness of the results 
obtained by the continuation procedure can be challenged immediately. 

Stewartson (ref. 15) studied the fluid motion induced by the impulsive motion of a 
semi-infinite plate (with velocity U 0 ) in its own plane (x > 0) and found that for the region 
0 = x = U Q t, two boundary conditions are needed at each end of the region because the sign 
of the convection term changes in this region. This sign change implies that small dis- 
turbances travel in both directions, x/U 0 t increasing and x/U 0 t decreasing. There- 
fore, the solution now depends on conditions at x = 0 and at x = U 0 t. The limiting 
behavior of the boundary -layer solution must be Blasius as x approaches 0 and Rayleigh 
as x approaches U 0 t. Lam (ref. 4) also arrived at the same conclusion in treating a 
more general problem of a shock moving past a semi-infinite plate (with the coordinate 
system fixed at the leading edge of the plate). He found that for this problem the Blasius 
profile must be used at x = 0 and the Mirels (ref. 7) solution must be specified at 
x = U 0 t, where U 0 is the velocity of the fluid in the free stream following the shock. In 
the limiting case of a weak shock, Stewartson’ s and Lam’s problems are identical. 
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The complete boundary -layer problem on the wall of a shock tube also has similar 
characteristics. Once again two boundary conditions are needed both at x = 0 (diaphragm 
location) and at x = U 0 t (the contact discontinuity) to complete the solution in the inter- 
action region (0 = x = U Q t). Gupta (ref. 12) and Ban and Kuerti (ref. 13) have used the 
Mirels (ref. 8) solutions at x = 0 and x = U Q t in treating this problem. 

References 4 to 19, with the exception of Stewartson, have considered the transient 
boundary layers formed by the passage of a shock wave. Detailed literature relevant to 
the non-shock-induced time -dependent boundary layers may be found in references 20 
and 21. 


PROBLEM FORMULATION 
Governing Equations 

The two-dimensional time -dependent compressible laminar boundary -layer equa- 
tions for a binary mixture have been taken as the flow-governing equations for the problem 
under consideration. These equations, including viscous, heat conducting, and diffusional 
effects are (see refs. 22 and 23): 

Continuity of mass: 


+ — (pu) + -|r(pv) = 0 


9t + 9x 


ay 


(i) 


Continuity of species: 


3Ci 3Ci 


ac 


pi + U + V ) = (pc. v’.) 

at ax ay/ av' 1 1/ 


Conservation of momentum: 


ay 


(2) 


/ 9u au .3u\ a / 8u\ 

p — + u — + v — = — u — 

ax ay/ 9y\ ay/ 

Conservation of energy: 


(3) 


p ih + u 8h + v »\ auy . aj 
\at ax ay/ lay/ ay 


(4) 


where v'-^ in equation (2) and q in equation (4) are defined, respectively, as 

* - C 1 . D 1 T 3 ^ T 

V 1 12 8y pcj 9y 


(5) 
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-q = k — + D 12 

3y 12 


/_ 3 In c. 


3 In Cr 


P h i 


+ hr 


Dl2\ c l 


^2 

3 In T \ 

c 2 

9 y I 


‘1 ay T 2 ay 
The gas mixture is assumed to be ideal and calorically perfect; that is, 
p = pRT = P L + p 2 = pT^jRj + c 2 R 2 ) 

h = hj + h 2 = Cjhj + c 2 h 2 = T(c lCp)1 + c 2 c p ^ 


(6) 

(7) 

(8) 


The rectilinear coordinate system chosen is fixed to the flat plate where x is the 
distance along the surface measured from the leading edge, and the dimension y is 
measured normal to the surface as shown in figure 2. 



Location of 
leading-edge 
particles at 
time t 


Figure 2.- Schematic diagram of flat-plate model 
showing the coordinate system used. 


The following assumptions have been incorporated in this study: 

(1) The flow external to the boundary layer is unaffected by the boundary -layer flow. 

(2) The contact surface in the inviscid flow is assumed to be thin; therefore, no dif- 
fusion is allowed across the interface. This assumption is justified for a thin interface, 
that is, for the flow situations where the thickness of the interface is small compared with 
the length of the plate. Also, the relaxation of the boundary layer on the plate surface in 
the region where the Blasius state exists is not substantially affected by this restriction 
of zero-thickness interface because the approach to the Blasius value is asymptotic. 

(3) For short -duration gas flows over solid boundaries having high values of con- 
ductivity and heat capacity per unit volume, the wall temperatures may be assumed to be 
constant. This assumption has been used consistently in this work. 
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Limitations Required for Conical Similarity 

The inviscid flow in a shock or expansion tube is a function of the conical coordi- 
nate x'/t' as shown in figure 3. Since no viscous-inviscid interation is considered in 
the present analysis, the introduction of a flat plate alined with the stream at length £ 
from the secondary diaphragm does not alter the conical similarity of the inviscid flow. 
However, solution of the viscous flow over the flat plate requires specification of bound- 
ary conditions at the leading edge. Also, a time -similar flow may be specified if the 
coordinate system which originates at the leading edge is chosen. In this coordinate sys- 
tem the similar flow regions can now be specified. In a shock tube problem (x,t coor- 
dinate system), the flow between the shock and the interface is of interest. In this case 
the flow is invariant at constant x/t up to the time the interface crosses the leading 
edge of the model. However, for an expansion tube the flow of interest lies between the 
interface and the expansion fan. With the arrival of the interface, a new boundary -layer 
flow develops which again is time -similar if time is measured from the instant of arrival 
of the interface at the leading edge (x*,t* plane) provided a constant boundary condition 
can be specified at the interface x*/t* = U Q . In the general case, however, this condi- 
tion is not possible since the interface traverses the shock-initiated flow in the x,t plane. 



Secondary 

diaphragm 
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Figure 3-- Schematic representation of the test-section flow- 
in distance-time plane. 







(See fig. 3.) Since no general time-similar solution is possible, two limiting flow situa- 
tions have been treated: 

(1) If the interface arrives at the leading edge closely behind the shock wave, the 
time lag At’ may be taken as vanishingly small. The interface may be assumed to be 
coincident with the leading -edge particle line. (See fig. 4.) For this limit to be desig- 
nated the MR (or Mirels) limit, the x,t plane may be conveniently used since the x,t 
and x*,t* planes now coincide. 



Figure 4.- Distance-time diagram for the MR limit. Time delay represented 
by the shaded region is assumed to be vanishingly small. 

(2) On the other hand, if the shock wave is far ahead of the interface, that is, At' 
approaches «°, then the interface advances into a boundary layer which has almost com- 
pletely relaxed to a steady-state (or Blasius) boundary layer. This limit, to be called 
the BL (or Blasius) limit, may conveniently be treated in the x*,t* plane (see fig. 5) 
after obtaining the fraction of the model length ( a = -p-V from the x,t plane solutions; 

over this length the boundary layer behind the progressing shock wave may be assumed 
to be Blasius (BL). The value of a required to establish the BL limit may be obtained 
from a separate analysis similar to that by Cohen and Trimpi (ref. 6), Gupta (ref. 12), 
and Lam and Crocco (ref. 16). These analyses give the value of a~ 0.3 for most of the 
flow conditions considered herein. It should be noted for this limit that with the arrival 
of the interface, the accelerating-gas Blasius boundary layer is disturbed and swept out 
by the test gas. A new Blasius boundary layer of the test gas now develops over a section 
of the model length near the leading edge. 
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diaphragm 


Figure 5.- Distance-time diagram for the BL limit. 

Equations in Terms of Crocco's Conical Coordinates 

Since the Crocco coordinates put the governing equations in a very compact form 
and have the distinct advantage over others in that the singular nature of the equations is 
more strongly evidenced, this system has been employed for the present problem. The 
other commonly used transformation for this type of problem is Stewartson's transforma- 
tion. (See refs. 15 and 21.) However, the equations are quite cumbersome in this coor- 
dinate system. Moreover, the sign of the coefficient of convection term, which is both 
positive and negative so that propagation takes place in directions of a both increasing 

and-decreasing-.-is-not-known u ntil the velocity field is evalu ated. This problem has been 

discussed in detail in reference 21. ~ 

The governing equations (2) to (4) are first transformed by using the Crocco system 
in which the independent variables are 
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and the dependent variables are 
c = c(x,u,t ) 

T = 1 ■— = T (x,U,t") / (10) 

h = h(x,u,t) 

The independent variables (eqs. (9)) can be further recast (ref. 4) in the conical coordinate 
system. In this system the dimensionless independent variables are 


0 = — 

U 0 

p U 2 t 

- MV O 


and the dimensionless dependent variables are 

c(a,0 ,r) = c ^ 






H(of,0,y) = 


h-. - 

Hxfa.S.y) = -i— ! i 

n W J 

The conically self -similar solutions in the Crocco system can now be obtained if the fol- 
lowing relations are specified for c, !j), H, and C: 


\Jay 


M =0 
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Equations (2) to (4) with the use of expressions (9) to (13) finally become (see 
appendix A for details) 

Continuity of species: 

2 + \ 

3/3\Np r 3/3 ] 3/3 


'N 


Le 


Np, 


- 1 


9c 

9/3 


+ J 


= C«(/3- «)|g- 

da 


( 14 ) 


Conservation of momentum: 


c p 2 JL£ + ££ 0 = a(/3 - a)(c - — 
9^2 2 \ 8 a 9a 


(15) 


Conservation of energy: 



8 

1 9H 

9/3 

Np r 9/3 


N 


Pr 


( H i + 1 )< 1 - r )||+J(H 1 + i)(i -D 


+ 0 


80 

9/3 


N 


Pr 



+ 1) (i - n 




CO - a)a 21 (16) 


For the special case of a single-gas boundary layer (c = l), the terms containing Nl 6 
and J are removed from the system. 

The results presented in reference 24 indicate only a slight difference between the 
constant Np r and N Le profiles and the variable Np r and N Le profiles. There- 
fore, for simplicity Np r and Nl 6 will be assumed to be constant in the analysis with 
the values of 0.7 and 1.4, respectively. With constant Nl 6 and Np r , equations (14) 
to (16) reduce, after some algebraic manipulation, to 

Continuity of species: 


,2 A + *2 « + /^Le . \ i* £ + * M j . c(0 - «) £ 


N 


Pr 


9/3^ 


9/3 


i N 


Pr 


9)3 9/3 


d/3 


da 


(17) 


Conservation of momentum: 


v 9^2 2 v ^ '\ da v 9a 


( 1 . 8 ) — 
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Conservation of energy: 




+ (i - r)(H, + i)j* 


= C(/3 - O') Q! 


9H 

3 a 


(19) 


where 


(Hi + l) = (H + 1). 


>,l / C - P,2 


2il 

: P,2, 


+ (1 - c) 


BOUNDARY CONDITIONS 

Definitions of Various Boundary -Layer Regions 

The boundary -layer region in the BL limit may be divided in three parts (as shown 
in figs. 6 and 7): 

(1) Shock region M: 

UaSA 
0 ^ 0 ^ 1 
y = 0 



This region was, treated as suggested in references 4 and 16. 

(2) Limited interaction region (LI) (between c** = 1 and a = 1): 

(a* = 1) = a %. 1 
0 ^ j3 S 1 
y = 0 

The limited interaction here implies that this region is part of the accelerating-gas inter- 
action region (0 = a = 1) which existed prior to the arrival of the interface. Therefore, 
this region is influenced by the leading edge of the flat plate, but not by the test -gas region 
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Figure 6. - Schematic representation of the flow field in the 
BL limit (physical plane). 


0 


I Region 

LI Region ! M Region 

HI 

n j 1 

Boundary layer 
consists of a mixture 

Boundary layer and free stream | 

of test gas and 

both consist of accelerating gas j 

accelerating gas, 


whereas free stream 

1 

consists of test gas 

1 

only 

1 

1 


V U o 


Figure J.- Schematic representation of the flow field in the 
BL limit (a*, p-plane ) . 

which follows. The relaxation of this region, which occurs in a simple shock -initiated 
flow (such as in a shock tube) has been treated in references 4, 6, and 12. The analysis 
of reference 12 is used to define the "essentially steady-state" region and a value of 
a ~ 0.3 is obtained for almost all the flow problems analyzed herein. The value a ~ 0.3 
should be a sufficient guideline for determining the essentially steady-state region. How- 
ever , thi s value of a may be obtained with greater accuracy for a given flow problem 
by following the approach of reference 12. ' r - 

(3) Interaction region I (close to the leading edge): 

0 £ a* % 1 
0^)35 1 
y* ^ 0 
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The interaction region is bounded by two different steady-state flows. At the downstream 
boundary ( a * = 1), the interface advances into a boundary layer of the accelerating gas 
which is assumed to have relaxed to the steady state since the BL limit treats only the 
situation where a* = 1 lies at a s 0.3. At the upstream boundary, a* = 0, a steady 
boundary layer in the test gas exists. 

In the MR limit where the time interval between a = 1 and a* = 1 line 
approaches zero, the boundary -layer region may be divided in two parts (see figs. 8 and 9): 

(1) Shock region M: 


lJaSA 
0 £ 0 ^ 1 
y g 0 



The treatment of this region is the same as in the BL limit. 

(2) Interaction region I: 

0 s A = 1 

0 5 0 S 1 
y 5 0 

It should be noted that the limited interaction region LI in this limit is eliminated. 

With the definitions of the two limiting situations for the flow field in o;,/3 -plane 
given, attention will now be centered on specifying the boundary conditions. 



Figure 8.- Schematic representation of the flow field in the 
MR limit (physical plane). 
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Figure 9>- Schematic representation of the flow field in the 
MR limit (a, (3-plane) . 


Specification of Boundary Conditions 

The objective here is to specify boundary conditions to solve equations (17) to (19) 
for c, <f>, and H for 0 S (3 = 1 for the appropriate a range in the two limits: 

BL limit: 0 ^ a* ^ 1 

MR limit: 0 ^ a = 1 

In both of these limits, (p and H must be given as functions of j3 near the leading 
edge (a — 0, a* — 0) and at the boundaries of the shock region (a = 1) and the limited 
interaction region (a* = 1). It may be noted that solutions for the BL and MR limits differ 
because of the different boundary conditions that are specified on <p and H at a* = 1 
or a = 1. 

The boundary conditions will be specified in the following steps: 

(1) Shock and limited interaction region boundary conditions: The shock and lim- 
ited interaction region boundary conditions, needed for MR and BL limits, respectively, 

can be deduced from the same relations; — First-t-he-shock-region-is_consid ered. The 

solution in this region was established by Mirels (ref. 8) and Lam (ref. 4) and is briefly 
given as follows: 

By letting 

<P(a,l 3) =Q(a) M(/3) 
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and taking 


Q(a) 


= 

|A - a 


(A > 1) 


it is noted that for a = 1, <p = M(|3). Also, <f> is real for 1 ^ a ^ A. This relationship 
also implies no shear stress definition for a > A (where there is no flow). Therefore, 
for 1 £ a ^ A, 


M(0) 


(20) 


It may be noted that the shear stress (p — °° as a — A which is expected because the 
boundary -layer thickness goes to zero near a = A. 

By assuming H = Hm(/ 3), it is noted that C = C(H) = c(Hm) = C(/3) for 1 s a ^ A. 
With the use of these assumptions and equation (20), equations (18) and (19) yield, 
respectively, 


g2j 

dp' 


d*M C A - 0 _ 
>2 2 A - 1 


M 


1 9 H M 


\ N Pr dp 


Uo_\ _9M 9H M L _J_ 
2 h w j ‘ dp dp l Np r 


= 0 


(21a) 


(21b) 


The boundary conditions on M(/3) and Hm(/ 3) are 


9M 

dp 




= 0 


p=o 


M \p=l = ° 


(22a) 


J 


H M (/3 = 0) = 0 
H M (/9 = 1) = H 0 = 




(22b) 


l w 


Equations (21a) and (21b) are the Mirels (ref. 8) similarity equations in terms of Crocco 
variables. To specify the shock region boundary condition on concentration, it should be 
noted that the boundary layer as well as the free stream consists of one gas (accelerating 
gas) in the region 1 % a ^ A. Therefore, 

e(or = 1,0) = 0 (23) 

for all values of p. 
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The limited interaction (LI) boundary conditions are recovered by setting a =-a* 
and A = 0 in equations (21) to (23). Here A = 0 does not imply the absence of flow, 
but merely that the momentum equation (21a) is to be solved for the Blasius value of the 
shear function. By setting A = 0 in equation (21a), it reduces to the form identical to 
that of Blasius equation (in Crocco form) that is to be developed in the next section. 


(2) Leading -edge boundary conditions: The boundary condition at a* = 0 (BL limit) 
and at a = 0 (MR limit) will be provided by the solution of the steady flat -plate govern- 
ing equations. At a = 0, <j> = G(/3); and, since any approach to the steady -state condition 


must be asymptotic, 


3alo!=0 


= 0. The momentum equation (18) for such a case reduces to 


o£& + Sf =o 

9j3 2 z 


(24a) 


Similarly, the assumption of a single gas boundary layer at a = 0 and 

H Uo = H o<ffl 

and 


® =0 


9 « Iq !=0 


simplify equation (19) to 

\Npr 9/3 2 h w J 9 /3 9 /3 \ N p r / 
The boundary conditions on G(/3) and H q(/ 3) are 


(25b) 

Under the assumption of an all test -gas boundary layer at a = 0, the following con- 
centration boundary condition maybe“specified-for-all-values-of — /3: — - — - — — — — — - — - — — 

c (a =0,/3) = 1 (26) 

The BL limit boundary conditions at a* = 0 can be obtained by replacing a by a* 
in equations (24) to (26). 


H g (/3 = 0) = 0 
H g (/3 = 1) = 0 



(24b) 


(25a) 
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(3) Wall boundary conditions: In the MR limit boundary conditions on <p and H 
for a zero-slip, impermeable wall and constant wall temperature are 


H(a,0) = 0 


(05a51) 


(27) 


Since the plate surface represents a solid boundary, the following condition may be speci- 
fied on c: 

|§(a,/3 = 0) = 0 (OSagl) (28) 

d P 

This condition means that there are no normal concentration gradients of either gas on 
the surface of the plate since the plate surface is neither a source nor a sink for either 
species. 

The wall boundary conditions for the BL limit may be obtained from equations (27) 
and (28) by replacing a by a?*. 


(4) Boundary-layer edge or free -stream condition: The free -stream boundary con- 
ditions on <jt> and H for the MR limit are 


0(a,/3 = 1) = 0 

Wno * h w 

H(a,/3= 1) = 


(0 g a ^ if 

> 

(OSaSl) 


(29) 


In equations (29), the free -stream temperature has been taken equal to the plate tempera- 
ture. Also, the free stream up to the contact surface will consist of all test gas. Hence, 

c(a,0 = 1) = 1 (0 ^ a s 1) (30) 

Once again, the BL limit boundary conditions may be obtained by using a* in place 
of a in equations (29) and (30). 

The values of the parameters A, U 0 ^/h w , T 0 /t w , and so forth required in obtain- 
ing the boundary condition and for solving the governing equations have been obtained from 
the inviscid shock-tube relations (ref. 25). These relations are given in appendix B. 


FLUID PROPERTIES 


The fluid properties needed in the boundary -layer equations are p, Cp, c y , p, 
k, hj, h 2 , h, Dj 2 > and Di T . The properties developed here are those dealing with 
binary mixtures of nonreacting gases. The binary mixtures considered are 
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(1) Nitrogen (test gas) and helium (the accelerating gas). 

(2) Nitrogen (test gas) and argon (the accelerating gas). 

The mixture properties are functions of the mass fraction of the individual species and 
properties of the pure species comprising the mixture. Extensive calculations for the 
required transport and thermodynamic properties have been made and these properties 
are discussed in appendix C. It should be mentioned that the binary collision integral 
tl 11* 

which appears in the various transport properties expressions depends on the 
choice of the molecular interaction potential. In this work this collision integral evalu- 
ated by Monchick's (ref. 2) exponential potential has been used. This potential is known 
to be the true qualitative form of the repulsive intermole cular potential at high tempera- 
tures. As pointed out by Monchick (ref. 2), it seems doubtful that a simple Lennard- 
Jones potential can be used over more than part of the temperature range of interest. A 
comparison between the various transport properties, calculated by using the exponential 
potential and the Lennard -Jones potential, is included in appendix C. 

COMPUTING PROCEDURE AND APPLICATIONS 


Numerical Solution of the Governing Equations 


The boundary -value problem formulated in this analysis requires the solution of 
three coupled nonlinear partial differential equations (eqs. (17), (18), and (19)). Since it 
is necessary to specify boundary conditions both at the leading edge of the plate and at 
the interface, the present problem was solved by using the Gauss -Seidel line -relaxation 
method, a technique which is used in elliptic problems. (See ref. 26.) 


To obtain the finite -difference equivalent of equations (17) to (19), the derivatives 
appearing in these equations are replaced by the following central difference 
approximations: 


9F _ 

"F T? 

m+l,n ' r m-l,n 

+ o(h 2 ) 

da 

2h 

9F 

F F 

m,n+l m,n-l 

+ o(h 2 ) 

80 

2h 

a 2 F 

myn+T — ^ r m,n 

+ F 

- - -m^n 

9/3 2 

h 2 



-\ 


(31) 


The finite -difference molecule is shown in figure 10. 
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Figure 10.- Finite -difference molecule. 


By using these finite -difference expressions, equations (17) to (19) may be written 
as follows for solution along lines of constant 0. The superscript p refers to the 
iteration number. 





c 


P 

m+l,n 



cP- 1 + c p 
m,n+l m,n-l 



%r\ 

N Le/ 


^m,n+l " ^m,n-l 


^m,n 


P-1 

m,n+l 



hNp r / p_i jP~l \ h ^Pr jP-1 ^m,n+l ^m,n-l 
4 N Le \ m,n+l ' 4 m,n «£ m ^ 


(32) 
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a 


! m( a m “ ^ 

n) C m ,n 

tiP . ttP 

n m-l,n + n m,n 

“m(“m - S») C m ; 

H p 

n m+l,n 

2h 

I m,n_ 

2h Im,n 

U 2 0 2 

u o r m,n 

+ 1 

nLp-1 iH P V^Le" 1 ! 

- r 1 

hw Im,n 

Np r h 2 

: Im,n V m > n+1 

m . n -V N Pr 

r 4h 2 





r 2hA m , n m >" 


P-l 

m,n+l 


- H: 


m ,n-l 


I 


m,n 


Im,n\ m ’ n N Pr h 2 j 


(34) 


In equation (34), I m)n is given by 

2 

T _ 2 N Le ‘ 1 1 - T 1 ^m,n, _ \ N Le " 1 1 - T 

m,n = h^ r h2y n W +1 m - n Vl] ‘ Np r r 


<P 


x 


41,2 *m,n ^ m ."+l 


c m,n-l)(^ , m,n+l " ^m,n-l) + 


N Le ~ 1 (1 - D 2 1 

Np r r 2 4h 2 


<P 


HL£/ C -C \2 (1 - T) 2 1 ^m,n , u 

2 ( m,n+l m,n-lj + 2 2h *2 l m,n+l nijn-lj^ni^ 


x 


2h x 


m,n 


4>: 


x-M 
^•m.n 


( J 


\ 1 - T 1 ^m.n, 

m,n+l " °m,n-lj ' r 2h A m 


- J. 


f (‘W+l - ^m,n-l) J 


i -r j_ 
r 2h 


(35) 


Equations (32) to (34) have been uncoupled since only c p , 0 P , or H p at points 
(m-l,n), (m,n), and (m+l,n) is unknown in each of these equations, respectively. In addi- 
tion, all nonlinearities have been removed by using values of the dependent variables 
from the previous iteration in the coefficients. Each of equations (32) to (34) forms a set 
of ^ simultaneous difference equations for a given value of n when m varies 
from 2 to 1/h. Since the coefficient matrix of each of these equations is of the tridiagonal 
form, t he Thomas algorithm discussed in reference 26 provides a simple algorithm by 

which the solution may be obtained. To start~the numeri'ca'Esolutionpthe-mass-f-raction 

is assumed to be unity throughout the flow field as if the problem to be analyzed contains 
the same gas in the accelerating section as in the test section. With c = 1, equation (32) 
is identically satisfied. Next, after assuming an initial distribution of <p, equation (33) 
is solved by the Thomas algorithm along lines of constant /3, beginning at the wall and 
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continuing until the outer boundary is reached. At this point the pth iteration has been 
completed. The iteration of the momentum equation is continued until 

Max | <(P - (jP * | < e ^ 

where e j is of order of 10”^. This value of p is referred to as P max . With this 
solution to equation (33) and still assuming c = 1 , equation (34) is solved in the same 
manner p max times. With the solutions to equations (33) and (34), equation (32) is now 
solved iteratively p„ ov times with appropriate boundary conditions. The resulting 

lYlclX 

c -distribution is further used to improve the solutions for 0 and H. This process is 
repeated until convergence to the desired accuracy 

(q,0) - (fcHa, 0) 

< ^i” 1 ( Q ! »°) 

is obtained. Here e 0 is a small number taken as 10"®. This criterion for conver- 
gence was selected since <£ w is a sensitive parameter characteristic of a solution. A 
similar criterion is used in reference 27. 

In the present calculations h = 0.01 was used. To improve the convergence rate, 
the initial distributions were obtained by solving equations (32) through (34) for a coarse 
mesh (h = 0.1) and then with linear interpolation these coarse solutions served as initial 
guesses in the iteration with fine mesh (h = 0.01). This procedure significantly reduced 
the number of iterations needed for convergence to the desired accuracy. The procedure 
employed in the computer code is given in the flow chart (fig. 11). 

Evaluation of Numerical Accuracy 

The accuracy of the present numerical technique was evaluated by comparing with 
cases given by Lam (ref. 28). By assuming c = 1 for a single-gas boundary layer and 
C = 1 for the viscosity law, equations (18) and (19) can easily be converted into the equa- 
tions treated in reference 28 for a shock tube. The cases of reference 28 were computed 
by using the present numerical approach and a comparison is made in table I. There is 
very good agreement between the two results. It is thought that the present results are 
more accurate because of the use of double precision in the computations which has the 
effect of keeping the round-off errors to a minimum. 




Boundary -Layer Parameters 


For a binary mixture, heat transfer at the wall may be obtained from equation (6): 

V “ -t(f) w + ( h w,l - | T fll < 3 ' 
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TABLE I.- COMPARISON OF SHEAR PARAMETER COMPUTED BY PRESENT 
NUMERICAL METHOD WITH THAT COMPUTED BY LAM (REF. 28) 

[a = 1.2 and C = lj 


0 

<£(a,|3) from Lam's results for — 

<p(a,i 3) from present results for 

- 

a = 0 

a = 0.3 

a = 0.6 

a = 0.8 

a = 1.0 

a = 0 

a = 0.3 

a = 0.6 

a = 0.8 

a = 1.0 

0 

0.330737 

0.345968 

0.517542 

0.729052 

1.145952 

0.331981 

0.340742 

0.518983 

0.736626 

1.145524 

.2 

.328361 

.341031 

.496062 

.697242 

1.095612 

.329971 

.335141 

.496927 

.704008 

1.095701 

.4 

.314174 

.321737 

.435211 

.606915 

.953981 

.315757 

.315419 

.434982 

.611645 

.954121 

.6 

.274307 

.277838 

.338874 

.462942 

.727993 

.275827 

.272653 

.338593 

.465303 

.728219 

.8 

.189752 

.190712 

.209048 

.266147 

.418522 

.190831 

.188173 

.208274 

.276861 

.418886 

.99 

.0171208 

.0171100 

.0178980 

.0177443 

.0266531 

.0181601 

.0176282 

.0177190 

.0187913 

.0274053 


By introducing Np r and noting that 
9F _ ( f b p w u o 9F 

9y ” Mw a P 


and 


=0 


9/3 


(for the present problem) 


w 


from equation (36), after some algebraic manipulation, 


Pw U o h w^w 

1 9H 

1 1 - F 

qw 

Np r 9/3 ' 

c w + (1 - 


J 

_ w 


is obtained where 


R 


w ,x 


Pw_ u o x 

M- 


= ay 


w 


(37) 


Equation (37) may now be used to obtain the following expression for the Stanton number: 

(38) 


N St = 


c p 1 

w \ 

1 u p h. 

1 

9H + 

1 -r J 

('I r w,x) N2 (t^- 'jl 

( /i w p w) Nz Mn 2 

_Np r 

9/3 

c + (i - c)r 


w 


In order to evaluate the recovery temperature T r appearing in equation (38), a 
recovery factor is defined as follows: 


T r - Too 

T s t ” T°o 


(39) 


where T s t is the stagnation temperature in the free stream (having a temperature T ^ 
behind the interface. For the sake of simplicity, the recovery factor is assumed to be 
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related to the Prandtl number by 


r = N 


1/2 

Pr 


( 40 ) 


In the formulation, Too = T w and 


+ 


U r 


st 1 °° 2 


( c p)n 2 


Therefore, from equations (39) and (40) 


^-l = Ni / 2 
T v 


U f 


■ w 


Pr 2( h w) N2 

Similarly, the expression for the local skin -friction coefficient may be obtained as 

4>V 


(41) 


c f = 21 




w 


( P w M w) N2 ]/( R W,x) N , 


(42) 


In the limit of an all N 2 boundary layer, c is identically equal to 1 and expres- 
sions (38) and (42) reduce to, respectively, 

0 


N St " 


w 


^( R w,x) N2 \ 

2(p 


(T 

T w 


^ - 1 


l [am 

Np r \ 8 /3/w 


(43) 


c f = 


w 


]j( R w,x) 


(44) 


N 2 


From equations (41), (43), and (44), one can also write 
N st .(9H/8/3), 


'w 


^ N ?/ 2 U ° 2 


(45) 


Pr ( h w) 


N 2 


Results in Terms of Physical Variables 

The obtained numerical results may furthefbe'Utilized'to-examine-t-he-result-ing- 
behavior of the solution in terms of the physical variables. From the definition of the 
transformed variables, 
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^w U o 3 ^wV /2 


<t> 


or 


8U ,f p w U o 3 ^wV /2 0 

3y \ x / M 


so that 


V (a b) = ( 5 \ 1/2 f ^ £ l£ Jl d| 

y( ’ } ip w U 0 M w j J 0 *(«,« 


Now 


M 


= C 


M,„P 


w'w 


Ri h 

p = Cp w p w 


Mi 

c + (1 - c) == — 
M 2 


where 1/p was obtained from equation (7), after rearrangement, and 

i C JZ£ 


h h wl Cw + (l_Cw ^ 

iL = ^(H + 1) 

C P C P,1 C + (1-C)^i 

c p,l 


Hence, from equations (47), (48), and (49), the dimensionless ordinate is 

y(o?,/3) Pr ( q?) ^ 


rj(a,0) = 


V Un 


1/2 y. 


^ = f 

n - 1 K,i J ° 


^1 p 0 


where 


0(«,4) 


P R («,S) = C(a-,£)<| c(a,4) + (1 - c(a,£)]=i 




(46) 


(47) 


(48) 


(49) 


(50) 


x 


C i 

:(a,0) + [l - c(o, 
-Eil y (H + 1 ) 

c(a,£) + fl - c(q!^)]-2^ 

L JC P,1 


(51) 


' • To evaluate the dimensionless ordinate rjla,—\, the integral appearing in equa- 

\ u o / 

tion (50) is evaluated at fixed a for various 0 values by using Simpson's rule. This 
procedure gives profiles of v(&, at various a values which may be interpreted as 

velocity profiles ^-( 77 ) for various values of a. 
u o 
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From expression (50), 


y _ rj__ 


14 


■ 0.91 
3 

) 

Jm r («,I)/0(«,^ 

]d£ 


( 52 ) 


may also be written. The temperature in the boundary layer may be obtained from 
equation (49). 


DISCUSSION OF RESULTS 

Calculations of the mass fraction of test gas (nitrogen), skin friction, and heat- 
transfer coefficients for transient laminar boundary layer of a binary mixture of perfect 
gases on a flat plate in the interaction region of an expansion tube are presented in fig- 
ures 12 to 15 for a wide range of shock Mach numbers. The gas considered in the test 
section is nitrogen. The gases considered in the expansion section are either helium, 
nitrogen, or argon. Helium has been most commonly used as the accelerating gas in 
expansion -tube tests. Nitrogen and argon are, however, considered (with the same free- 
stream velocity U Q ) to determine the effect of molecular weight of the accelerating gas 
on the relaxation of the boundary layer. It should be noted that M s = 8.55 in argon, 

M s = 7.17 in nitrogen, and M s = 3 in helium give the same free-stream velocity of 
2.04 km/sec. In other cases considered, M s = 20.8 in nitrogen and M s = 8 in helium, 
both yield a free-stream velocity of 6.03 km/sec, whereas M s = 26.1 in nitrogen and 
M s = 10 in helium result in a velocity of 7.56 km/sec. Although the results for the 
higher Mach numbers are academic because of the assumptions made concerning the vari- 
ation of properties of nitrogen at high temperatures, these calculations were performed 
in order to indicate trends. Thus, although nitrogen is treated as a perfect gas in the 
binary mixture, its viscosity is obtained from the real-gas values of Yos. (See ref. 29.) 

In the present calculations, the plate surface was assumed to be at a constant tem- 
perature of 300° K and for convenience the free stream behind the interface was assumed 
to be at the same temperature. This last assumption is relevant for expansion -tube flow 
duplication but is inappropriate for simulation studies. (See ref. 30.) Perfect -gas rela- 
tions, which are given in appendix B, were used to obtain the conditions behind the shock 
for a given shock Mach number in the accelerating gas. 

Fig ure 12 sh ows distri butions of the mass fraction of nitrogen between the flat -plate 

leading edge and the interface for various values of shock Mach number in different gases" 
in the expansion section of the expansion tube for both the BL and MR limits. From these 
figures it may be concluded that the helium -nitrogen boundary layer has relaxed to the test 
gas (nitrogen) for values of a up to 0.3. For example, it relaxes to about 95.5 percent 
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nitrogen at the wall for M s = 3 in helium, to 96.5 percent nitrogen at the wall for 
M s = 8 in helium (with thermal diffusion), to 97 percent nitrogen at the wall for M s = 10 
in helium (with thermal diffusion); and to 98 percent nitrogen at the wall for M s = 8.55 
in argon. For the BL limit, these percentages are about 1 to 2 percent lower. These 
figures also indicate a physically plausible trend that the effect of thermal diffusion 
increases with the increase of shock Mach number. Figure 13 shows the relaxation of 
the binary boundary layer to a nitrogen boundary layer as the shock Mach number is 
increased. (The data points in this figure have been connected by straight lines.) It is 
revealed that the binary -gas boundary layer relaxes faster on the plate surface (/3 = 0) 
in terms of the parameter a (or a;*) as the shock Mach number in helium is increased 
from M s = 3 to M s = 10. This trend is consistent for the MR as well as the BL limits. 


Figures 14(a) and 14(b) give the Stanton number distribution for the MR limit. Fig- 
ure 14(a) also includes one curve (for M s = 3 in helium) which gives the Stanton number 
^still multiplied by (l/Rw^^j f° r helium considered as the test gas. This figure indi- 
cates that for a ~ 0.3, the Stanton number relaxes to within 6 percent of the steady-state 
nitrogen boundary -layer value. In figure 14(b) Stanton number distributions are com- 
pared for nitrogen (with M s = 26.1) and helium (with M s = 10) which have the same 
free -stream velocity U Q . This figure also indicates that for a ~0.3, the Stanton num- 
ber has relaxed to the steady -state nitrogen value (within about 9 percent). Also included 
in this figure is the curve for M s = 26.1 in nitrogen when C (defining the viscosity 
law) is assumed to be unity. It may be noted from these solutions that when the gas used 
in the expansion section is a light gas (such as helium), the boundary layer relaxes faster 
to the steady-state test -gas condition. However, a different trend is indicated for the BL 
limit in figures 14(c) and 14(d). In this case a light accelerating gas results in a slower 
relaxation when compared with the other gases. A dip is noticeable in figure 14 for 


With thermal diffusion 


Limits 




MR ) 
MR l 
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BL J 


c * 0.95 
& = 0 


c = 0.99 
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Figure 13 •- Effect of shock Mach number on relaxation rate 
for nitrogen-helium boundary layer . 
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Figure l4.- Stanton number distribution. 
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Figure lh.- Continued. 








cases when a lighter gas is used as the accelerating gas. The effect of thermal diffusion 
is less noticeable in the Stanton number distribution as compared with the nitrogen mass 
fraction distribution. 

In figure 15, distributions of the local skin -friction coefficient are given for the MR 
and BL limits. In figures 15(a) and 15(b), which are for the MR limit, an interesting 
result is noted. These figures show that when helium is the accelerating gas, the skin- 
friction coefficient deviates only slightly from the steady-state value. This result may 
have significance in experiments since a faster instrumentation response to the steady- 
state value may be expected. For the other cases, a ~ 0.3 is the location where the 
skin-friction coefficient has relaxed to within 5 percent of the steady-state skin-friction 
value in the test gas. The skin -friction coefficient for the case when helium is consid- 
ered as the test gas (for Ms = 3 in helium) is also depicted in figure 15(a). 

The skin -friction coefficient for the BL limit is shown in figures 15(c) and 15(d). 

In the BL limit, however, using the test gas as the accelerating gas appears to give faster 
relaxation. Once again, a* = 0.3 is the location on the plate surface upstream of which 
the skin-friction coefficient has relaxed to the test-gas steady-state value. 

A qualitative comparison with the integral solutions of Trimpi and Cohen (ref. 6) 
for a shock tube using two different fluids can now be made. Their results show that for 
a < a' crit , most of the boundary layer is comprised of driver gas. This value of Q! cr it 
corresponds to a = 0.3 obtained here. For a five- and six-term velocity profile, they 
found the value of 0.4 for af cr jj-. The comparison is approximate because the integral 
results do not give the exact extent of the driver -gas boundary layer. 

Figures 16 to 18 give the distribution of velocity, temperature, and mass fraction 
of nitrogen, respectively, through the boundary layer as a function of the physical vari- 
able y/6. The results presented in these figures are for M s = 8 in helium in the MR 
limit. The velocity profiles are smooth and validate the assumption implicit in the 
Crocco transformation that u(x,y,t) is a monotonically increasing function of y through 
the boundary layer. The velocity profiles at a = 0 and a = 1 are the self -similar 
Blasius and Mirels (ref. 8) profiles, respectively. 

It should be noted here that 6 represents the velocity boundary -layer thickness, 
evaluated by assuming that the edge of the velocity boundary layer is at /3 = 0.99 (that is, 
at the point where the velocity is 99 percent of the free -stream value). The thermal and 
the concentration boundary -layer thicknesses will be much larger and may be evaluated 
by using the value of 0 much closer to 1 (for example. fi = 0 . 9.9.9).— — T 

Finally, figure 19 gives the relaxation of transformed boundary -layer thickness 
from Mirels' value at a = 1 in helium (for M s = 8) to the Blasius value at a = 0 in 
nitrogen. Once again, only the MR limit result is shown. A rapid steepening of the 
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Ng-He boundary layer with 
thermal diffusion 



(a) M s = 3, M s = 8 in helium; M s = 7*17 > M s = 20.8 in nitrogen; 
M s = 8.55 in argon; MR limit. 

Figure 15 . - Distribution of local skin-friction coefficient. 
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■ Ng-He boundary layer with 
thermal diffusion 

Ng-He boundary layer without / 

thermal diffusion / 

Ng-Ng boundary layer 


/ 


/ 


/ V_ M = 26.1 in N 2 (c = 1 ) 



I I I I I I 1 I I 
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Nondimensional distance-time variable, a 


10 in helium; M s = 26.1 in nitrogen; MR limit. 
Figure 15-- Continued. 



N 2 -He boundary layer with 

thermal diffusion 

N 2 -He and Ng-Ar boundary layer 

without thermal diffusion 

n 2 -n 2 boundar y la y er 



= M s = 8 in helium; M s = 7-17, M s = 20.8 in nitrogen 
Ms = 8.55 in argon; BL limit. 


Ng-He boundary layer with 
thermal diffusion 
Ng-He boundary layer without 
thermal diffusion 
N 2 -N 2 boundary layer 



M s = 10 in helium; M s = 26.1 in nitrogen; BL limit. 
Figure 15.- Concluded. 






Figure 18.- Mass fraction of nitrogen through the boundary layer 
for M s = 8 in helium. MR limit. 



Figure 19.- Relaxation of boundary-layer thickness to Blasuis value 


(Mo, 


for M s = 8 in helium with nitrogen as the test gas. MR limit. 



curve may be noted for a — 1. This result is due to the temperature jump contained in 
the thermal boundary condition at the singular point a = j. 3 = 1.0. As expected, the 
boundary -layer thickness approaches zero (and the shear function approaches infinity) at 
the foot of the shock which in this case is at a = A = 1.356. 

CONCLUDING REMARKS 

A theoretical investigation of the relaxation of accelerating gas to the test gas in the 
boundary layer over a flat plate mounted in an expansion tube has been conducted. The 
time -dependent compressible laminar boundary -layer equations of a perfect binary gas 
mixture have been taken as the flow -governing equations. Complete mass and thermal dif- 
fusion effects for the binary mixture in the boundary layer have been considered. Self- 
similar solutions of the governing equations expressed in Crocco variables have been 
obtained by the Gauss -Seidel line -relaxation method. The results obtained correspond to 
two limiting cases in which the flow external to the boundary layer on a flat plate can be 
considered to be similar in time or "conical." These conically similar limits are (1) when 
the time lag between the arrival of the shock and the interface at the leading edge of the 
plate is very large; and (2) when this time lag is vanishingly small. 

The following assumptions have been made in the present analysis: gases have been 
treated as the binary mixture of perfect gases, the real -gas viscosity has been used for 
nitrogen, and the interface has been assumed to be of zero thickness. This last assump- 
tion is justified for a thin interface, that is, for the flow situations where the thickness of 
the interface is small compared with the length of the plate. Also, the relaxation of the 
boundary layer on the plate surface in the region where the Blasius state exists is not sub- 
stantially affected by this restriction of zero -thickness interface because the approach to 
the Blasius value is asymptotic. 

The fluid properties required in the numerical evaluation of the equations have been 
obtained in detail from the exponential potential function, which is considered to be more 
accurate than'the Lennard-Jones potential at high temperatures. A comparison between 
the properties calculated from these two potential functions at high temperatures is 
presented. 

The numerical results obtained are for shock Mach numbers ranging from 3 to 10 
in helium (considered as the accelerating gas). The gas considered in the test section is 
nitrogen. ■ Argon and nitrogen are also considered as the accelerating gases for shock 
— Mach-numbe rs~giving equivalent velocities iiTthe'free 'stream behindThe shock. The 
results obtained predict the time required for an all-helium or an all-argon boundary layer 
to relax to an all-nitrogen boundary layer at any point on the plate. The results indicate 
that a steady-state boundary layer containing more than 95-percent test gas ("perfect" 
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nitrogen) exists over a plate length equal to about three -tenths of the distance traveled 
by the interface in the free stream from the leading edge. The Stanton number and the 
local skin-friction coefficient, however, relax to within about 10 percent of the steady- 
state test -gas value for about four -tenths of the plate surface. 

Using a lighter gas (such as helium) as the accelerating gas increases the test 
time for the case when the shock wave and the interface arrive on the model almost 
together. However, using a heavier gas (such as argon or nitrogen) as the accelerating 
gas increases the test time for the case when the interface lags behind the shock wave 
substantially. 

Finally, since the conversion of the boundary layer is near completion for about 
three-tenths of the plate length in the two limiting situations analyzed here, it is thought 
that all intermediate cases should lie within the range of the results presented in this 
report. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 4, 1972. 
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APPENDIX A 


TRANSFORMATION OF FLOW-GOVERNING EQUATIONS 


The Crocco dependent variables are 
c = c(x,u,t) 

T = M r(x,u,t)> 

h = h(x,u,t ) 


and the independent variables are 
x s x 

u s u(x,y,t) 
t = t 


(Al) 


(A 2) 


For a function F(x,y,t) transformed to F(x,u,t), the transformation equations for its 
derivatives will be 

9F _ _9F + _9F 9u q 

9x 9x 9 u 9x 


9F _ 9F 9u 
9y 9u 9y 


(A4) 


9F _ 9F 9u 9F 
9t ” 9u 9t + 8t 


(A 5) 


If F = y, then 

.& = 0 ^ 
9x 

= i 

ay 

2L= o 
at 


(A6) 


Er.om_equations- (A3)_and-(A6)-,-it-follows-that 
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or 

9u _ 9y/ 9x _ y x 

" " 9y/ 9u “ ‘ y a 

Similarly, from equations (A5) and (A6) 

is. = 0 = + 

9t 3u 91 at 

or 

9u y f 

ot--y. 

Also, from equations (A4) and (A6) 

1 = 

Qy au 9y 


(A7) 


(A8) 


or 

9u _ 1 _ T_ 

0y 9y/ 8u M 


(A9) 


By using equations (A7), (A8), and (A9) in equations (A3), (A4), and (A5), 


9F _ 

9F 

y x 

aF 

9x 

9x 

y- 

u 

9u 

9F _ 

r 9F 



9y 

4 9u 



9F 

9F 

y t 

9F 

at ‘ 

at" 

y U 

aa 


(A10) 

(All) 

(A12) 


In terms of Crocco's variables, equations (2), (3), and (4) may now be written as 


9c - 9c \_ 2 (9 ^Le 9c 


9t 


9x 


+ ■— -r + 


9u Np r 9u 



+ i 

pL( T 2) 


/ N Le A 9c , d 1 T 1 9T 

2 

|9u 


yNpr J 9u T M &u 


(A13) 



+ u 



(A 14) 
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APPENDIX A - Continued 



N Le ~ 1 3c + d 1 T 1_ 3t\ 
V N Pr 9u + T ^ 9u/ 


(A15) 


In obtaining equations (A13), (A14), and (A15), the overall continuity equation (1) has been 
utilized. 

Next, these equations are transformed to the conical coordinate system. In this 
system the dimensionless independent variables are 



This transformation, although completely general, is subsequently applied to a restricted 
class of flows in which the free -stream quantities are functions of a alone in order to 
reduce the complexity of the equations while still retaining all the features necessary for 
applications to shock-generated flows. 

The transformation equations for the derivatives are 

9F _ J_ _9F 

9x U 0 t 9Q! 

9F _ 1 9F 

9u U 0 8 /3 

tt 2 

— 9F— - — Ci!-9F Wq — 

at t da 
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APPENDIX A - Continued 

Equations (A13), (A14), and (A15) may now be written, by using equations (A17), as 

^Le 


T 2 3 


Le ^ + J 


Uq 2 9/3\N pr 3£ 


2U n ^L‘ 


— (r 2 ) 
3/3 ' ' 


\ N Pi 


'¥ (J 


u, 


= pp 


o_ 3c _ a - /3 3c 
^W 3y f 3a i 


(A18) 


t 2 3 2 t , , (ot-B 3t U 0 2 a T \ _ a 
U Q 2 3/3 2 \ t 9a ' "w 




3a 


^w 9y 


(A19) 


U 2 1 0 9 /3 

u o 1 


3 

1 3h 

9/3 

Np r 3/3 


N t 


s Pr 


||h 1 (i-r) + jh 1 (i-r) 


u . 

_ -(r 2 )" 

/ W 0 2 

3/3 ' T ' 


X 


N 


1 - h l( i - r)^ Le ~ + j) 


Pr 


Np r 9£ 


(A20) 

\^w 9y t da 


The following dimensionless dependent variables are now introduced: 
c(«,/3 ,y) = c 
M°t,p,y) = 


P w^ 


H(a,/3,y) = 


h - h 


w 


■w 


H 1 (»A7) = b_' hw 


‘w 




With the introduction of equations (A21), equations (A18) to (A20) become 


9/3 \Np r 3/3 ] '3/3 


Li. 1 |£ + J 


■Jpr 


9/3 


_ c /3c _ a - /3 3c N 
Sf ” y 9a 


(A21) 


(A22) 


-f i± + cf^-2 ii - ^ = 5/2^1 jc . ^ 

3/32 \ y 9a 3yy \ y 9a dy J 


(A23) 
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-2 JU r 
0 


*w 


a 

1 9H 

9/3 

Np r 9/3 


N Le ~ 1 
Npi 


( H 1 + i)(i -r)|| + j(H 1 + i)(i -r) 




+ ( Hl + r>(^ 


Le_1 ^ + J 


N 


Pr 


9(3 


= C|— - a ' P — 

i 9y y da 


(A24) 


The following form for <fi is introduced so that the momentum equation (A23) 
becomes independent of y: 

(A25) 

v«y v^-w ,x 

For energy equation (A24) , the following relationship is assumed: 

H(cK,/3,y) = R(y) N(a,/3) 

Now R(y) is specified so that the lim H is bounded and equation (A24) is independent 

y— 0 

of y. Such a requirement for R(y) means that R(y) = Constant and without any loss 
of generality this constant may be taken as unity. 

Hence, 9H/9y = 0 and implies that 

H = H(a,/3) (A 2 6) 

Equation (A26) is the result given by Lam (ref. 4). ^It may also be noted from eq. (A25) 

that lim <p = 0.) By following an argument similar to that for H, it is found that 
y— 0 J 

c = c(Q!,/3) 

8c .o r (A27) 

fy ~ J 

Also, since C = C(c,H), 


9C = 0 

dy 


(A28) 


By utilizing expressions (A25) to (A28) in equations (A22) to (A24), the following 
form of the governing equations is obtained: 


< P 


2 9 / N Le 9c T \ , 90 


Le i)£ + j 


;8ipr — j-w 




_(A29.)_ 


8(3^ * 


da 


da 


(A 30) 
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_d_ 

9/3 


1 3H 
Np r 9/3 


N 


- 1 


^p 1 + i)a-n|+j( Hl 




+ 


i)(i - r) 



9c 

9/3 



CO - a)a 


(A31) 



APPENDIX B 


SHOCK PARAMETERS IN AN IDEAL GAS 


To be consistent with the rest of the analysis and for the sake of simplicity, the 
following properties of the gas are assumed: 


(1) It obeys the ideal equation of state. 

(2) y - Cp/c v = Constant, cp = Constant. 

With these assumptions, the following relations are obtained from reference 25: 


A - — - y + - ^s 2 

Uo 2 M g 2 - 1 

U ° 2 _ 4 (y - 1 ) T 1 /„ _ _ 1_\ 2 
h o (y + 1)2 1 0 \ S 


(Bl) 


(B2) 


To 

T 1 


1 + 


2 (y - 1) 
(y + l) 2 


yM s 2 


J— + (1 - y) 
M s 2 


(B3) 


where T-^ is the temperature of gas in front of the shock. The temperature Tj has 
been taken as T w (=300 K) in this analysis. 


Further, 



Up 2 h 0 
h o ^w 




(B4) 


since Cp is assumed to be constant and T^ = T w . 
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APPENDIX C 


THERMODYNAMIC AND TRANSPORT PROPERTIES OF 
N2-He AND N 2 -Ar MIXTURES 


Thermodynamic Properties 

N 2 -He mixture.- In evaluating the density p of the binary mixture, the individual 
species comprising the gas mixture are assumed to be calorically perfect; that is, the 
specific heats of the pure species are constant. The density of the mixture of perfect 
gases is obtained from equation (7): 


T £ c,Ri 

i-1,2 

or 


p = ^~ 

R-, T 


(Cl) 


1 A Mi 

c + (1 - C ) 

M 2 

where Rj is the gas constant for nitrogen and MjyMg = 7 for the N 2 -He binary 
mixture. 

The specific heats for the mixture are given by 
c v = c.c Vfl + (1 - c)c v 2 

c P = cc p,l + ( 1 c ^ c p,2 

For calorically perfect He and N 2 , the following relation is also employed: 

c Zr 

c p,l _ 2 K 1 _ _l-4_ _ Qj2 
c p,2 |R 2 Mj/m 2 

Values of Cp 2 , c y 2’ anc * ^2 ^ or helium are 

c n = 5.1988 x 10 3 • m „ 2 
p ’ 2 sec2-K 

c o = 3.1189 x 10 3 - 
v ’ 2 sec2-K 


R 2 = 2.0799 x 10 3 


m‘ 


sec2-K 


(C2) 

(C3) 


(C4) 


(C5) 
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N 2 -Ar mixture.- Equations (Cl) to (C3) may also be used for obtaining p, c p , 

and c v for the N 2 -Ar mixture. However, c will represent the mass fraction of Ar, 
and Rj will be the gas constant for argon. Other quantities of interest for this case 
will be 




^i= 1.429 
M 2 

5n 5 R’ 

i_L = 2_40 = 0 c 

C P,2'7r "IK! 

2 z 2 28 


c n = 5.2032 X 10 2 — ^ — 
P’ 1 sec 2 -K 


(C6) 


= 3.1171 x 10 


2 m" 


v,l sec 2 -K 

m2 


Rj = 2.0819 X 10' 


sec 


2-K 




Here, R’ is the universal gas constant. 

Transport Properties 

The viscosity of the mixture is given by Wilke's approximate formula (ref. 31) 
applied to a binary system 


4' 


M = 


1 


P' 


x 2 x x 

1 + G 12 1+G 21^ 


(C7) 


where p-^ and p 2 are the viscosities of the individual species forming the binary mix- 
ture and 


G ij = 


'"ii l M, 


,^L Mi) 

’ v2 ll-+-=- | 
\ M 


TT\l/2 
1. 

V 


(i, 3 = 1,2) 


(C8) 
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x l = 


x 2 ” 


M 1 

c + (1 - c) =r - 
M 2 

1 - c 

+ (1 ■ c) 

m 1 /m 2 


N 2 -He mixture.- The viscosities of N2 and He were obtained as follows: 
For 300 K ^ T = 2000 K (based on ref. 32): 

ix N 0 = 5.7925 X i()- 7 T 0 * 613 N ~^ £ 

1 m 2 

For 300 K ^ T = 1300 K (based on ref. 33): 

p He = 5.0236 X 10 -7 t °- 64 7 N-S£C 

m 2 

For 2000 K = T § 10 000 K (based on ref. 34): 


(C9) 


(ClOa) 


(ClOb) 


M n = 1.2802 X 10~ 7 T°- 811 £:g|? 

2 m 2 

For 1300 K £ T S 10 000 K (based on ref. 34): 


(Cl la) 


M 


He 


1.2858 x 10 V - 839 N " sec 


rm 


(Cllb) 


To simplify the present analysis, a mean viscosity -temperature relationship was 
applied to both gases. This relationship was obtained with the help of equations (CIO) 
and (Cll) and references 32 to 34 as shown in figure 20. The mean viscosity formula 
obtained is 


t'mean= 5 - 4085 xl0 ' 7T °' 63?i ^f 


Mmean = 1 - 2908 x 10" 7 T°- 825 N ' sec 

m 2 


■x 


(300 K IT ^ 1550 K) 


(1550 K s t-s 10 000 K) 


} (C12) 


The viscosity obtained from equations (Cl 2) has been fitted to polynomials in tem- 
perature by using the method of least squares 




mean 


= A + BT + CT 2 + DT 3 + ET 4 + FT 5 


(300 K ^ T = 10 000 K) (C13a) 


m‘ 
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Temperature, °K 


Figure 20.- Temperature variation of helium and nitrogen viscosities and mean value used. 


where 

~N 

A = 1.29908164 X 10~ 5 
B = 2.83599596 X 10' 8 
C = -1.63078339 X 10" 13 l 
D = -7.03066415 x 10" 17 f 
E = 6.29878349 x 10“ 21 
F = -1.53160089 X 10 -25 


(Cl 3b) 


This fitting has been done to remove the discontinuity in the slope of the viscosity curve 
at 1550 K. The M me an obtained from these polynomials is displayed in figure 21. 

For temperatures higher than 10 000 K, the viscosity of N£ may be obtained from 
the work of Yos (ref. 29). These data have been fitted to polynomials in temperature: 

M = A + BT + CT 2 + DT 3 + ET 4 + FT 5 + GT 6 -- ~ g | - (C14a) 

2 m 2 
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Figure 21.- Polynomial fits in temperature to various viscosities. 


where 

For 300 K = T ^ 10 000 K: 

A = 3.01575924 X 10~ 5 
B = -1.3942708 x 10~ 8 
C = 2.94077531 X 10 -11 
D = -1.00839516 x 10 -14 ) 
E = 1.6046776 x 10-* 8 
F = -1.15742765 x 10~ 22 
G = 3.03739295 x 10’ 27 

For 10 000 K T = 20 000 K: 

A = 4.28597628 X 10 -2 
B = -1.85480431 X 10' 5 
C = 3.27841650 X 10' 9 
D = -3.00995582 X 10" 13 > 
. E = 1.51565354 x 10' 17 
F = -3.98056198 x 10" 22 
G = 4.27318273 x 10" 27 


(C14b) 


(Cl 4c) 
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For 20 000 K ^ T ^ 30 000 K: 

A = -3.23611473 x 10“ 2 
B = 7.97818283 x 10“ 6 
C = -8.06518064 x 10' 10 

D = 4.27745175 X 10“ 14 > (C14d) 

E = -1.25286254 x 10" 18 
F = 1.91811505 X 10“ 23 
G = -1.19731314 x 10 -28 

y 

Figure 21 shows the curve for viscosity obtained from these polynomials. This figure 
also contains the extrapolation of this polynomial fit to 42 000 K. This extrapolation was 
used in the evaluation of boundary -layer parameters for shock Mach number of 26.1 in 
nitrogen, since no viscosity data exist for nitrogen for temperatures higher than 30 000 K. 

The thermal conductivities of the individual species are given by the Eucken rela- 
tion (ref. 22, p. 499) 

k i = lfe - 5 ) c v,i^i (i = N 2 , He) (C15) 


The mixture conductivity, as shown in reference 35, is obtained by replacing the viscosi- 
ties of the species, / and p 2 , in equation (C7) with the conductivities k^ and k 2 
given by equation (Cl 5). The quantities given by equation (C8) are still evaluated 

by using the pure species viscosities (ref. 35). 


The binary diffusion coefficient is computed by the expression (see ref. 22, p. 539) 


h3/2 


pD 12 = 0.002628 


~(mi + m 2 ) 


2M 1 M 2 


1/2 


atm -cm 2 


2 

CT 12 n i2 




sec 


(C16) 


where a^ 2 is some length parameter defined by the potential function 0. 


(1 11 * 

The collision integral fi^ 2 ’ ' in expression (C16) depends on the choice of molec- 
ular interaction potential. This integral has been evaluated by Monchick (ref. 2) on the 
basis of the potential 

0 = Ae" r//p (C17) 


where A and p are (ref. 34) 



A, J 

p, m 

He - He 

6.18411 x 10' 17 

2.20 X 10 -11 

n 2 - n 2 

2.16284 x 10' 15 

2.63 X 10' 11 
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For the mixtures, the following combination rules (ref. 34) are used: 


a i 2 - (A^y 




p 12 2 \ P 1 p 2 


The parameters which are most useful in fitting transport properties are p 
and a, where 

o = In A. (C19) 

Here k is the Boltzmann gas constant and T is the absolute temperature in kelvins. 
By following Monchick (ref. 2), the integrals of interest to gas transport theory may be 
put in the form 

n (l,i) = /j£T\ 1//2 J dx x 2 e _x Q^^(x) (C20) 


Q^(x) = 27T ^ j^l - cos x(x)jb db 


x = 

2kT 


where p is the reduced mass, g the relative velocity, and b the impact parameter. 
For a given 0 (intermolecular potential), x is given by 


= ’- 2 I 


r -2 b dr 


r ° (i - _ 20 \*/ 2 

\ r 2 Mg 2 / 

Here r Q is the distance of closest approach. 

With certain simplifying assumptions (ref. 2), x ma y be rewritten as 
X = 7r - 4k’ e Q 


(C23) 


(C24) 


k’ = 


cos 2 e 0 

cos 2 0 


sin e c 
sin 6 


(C25) 
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. . -1 b 

6 0 = sin — 


0 = sin 1 £ 


(C26) 




Z = - 


a - In x 


u = cos e r 


z 2 = 1 - ± 


(C27) 


With the help of equation (C24), equation (C21) may be simplified to 

Q <i) = §ie!i’ l(x) 


(C28) 


where 


I^x) = du(l - cos x) <jy + 4 -j - u 3 (4 In u + 1) 


+ V 


4 In u(2u 3 lnu+u 3 - J 


(C29) 


By utilizing equation (C28), expression (C20) may finally be put in the form 

a/2 




where 


and 


\ 1 t 3 (-4 In t) 2 

^ ln(-4 In t) 

J 0 

a 


I\ dt 


(C30) 


(C31) 


t = e-/ 4 

Monchick (ref. 2) has provided the tabulation of integral 1^ ^ as a function of a. 

Hirschfelder et al. (ref. 22, p. 526) have given their formulas for transport proper- 
ties for the first approximation in terms of the quantities 


a 2 0 (M>* - 
a 12“l2 


/'2ttp') 1 / 2 j 2 (1,1) 


UT 


f (2) - 


(C32) 
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Substituting for fi' 


( 1 , 1 ) 


from equation (C30), equation (C32) yields 


~2 r»(l.l)* . ^Al.l) 

a 12“l2 21 


(C33) 


By using equation (C33) in expression (C16), the following relation is obtained for the 
binary diffusion coefficient: 


4.3822 X 10' 2 T 3 / 2 „ , 
pD 12 5 - N/sec 


a 2 1 


( 1 , 1 ) 


(C34) 


where T is the absolute temperature in kelvins. Figure 22 shows the plot of pDi 2 as 
a function of temperature. Polynomial representation of this function is given as follows: 

PH 12 = A + BT + CT 2 + DT 3 + ET 4 + FT 5 N/sec (300 K ^ T ^ 10 000 K) (C35a) 


where 

A = -4.20571429 x 10 ° ^ 
B = 2.48592499 x 10 -2 
C = 4.10466471 x 10 ' 5 I 
D = -2.16644279 x 10 -9 ( 
E = 1.49001648 x 10 “ 13 
F = -4.43803342 X lO' 18 ^ 


(C35b) 


If, instead of using the Monchick’s potential to obtain and pD , the 

Lennard -Jones ( 6 - 12 ) potential is employed, the coefficients for the polynomial 2 fit 
(eq. (C35a)) become: 


A = -5.13254756 X 10° 

B = 2.95342073 x 10 " 2 
C = 3.14445935 X 10 ‘ 5 
D = -1.93226973 X 10 ' 9 / 
E = 1.27828289 X 10 “ 13 
F = -3.92041533 X 10 ' 13 


(C35c) 


For Lennard-Jones (6-12) potential, the following force constants are required 
(ref. 22 ): 



a, m 

e/K,K 

He - He 

2.576 X 10 ‘ 10 

10.22 

n 2 -n 2 

3.749 X 10 " 10 

79.80 


where e//c is the molecular potential energy parameter. 
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For the binary mixture, the following combination laws (ref. 22) are employed: 


CT 12 “ 


a l + a 2 


e 12 = V € " € 


12 


(C36) 


where the subscript 

(1,D* 


n. 


12 


1 refers to a heavy gas, subscript 2 to a light gas. The integrals 
for the Lennard -Jones potential are tabulated in reference 22 as functions of 


reduced temperature T* (= kT/c^- Although these tabulations give for tem- 

peratures greater than 10 000 K, pD^ obtained from equation (C16) may not be very 
accurate beyond 1000 K, the reason being that the force constants required in equa- 


tion~(C16)~are 'bbtaihecTfrom viscosity~data generally in the range from 300 K to 1000 K. 
For the sake of comparison the values of pD^ obtained from Monchick's exponential 
potential and Lennard-Jones (6-12) potential have been included in figure 22. The two 
molecular interaction potentials seem to give the same values for pD^ f° r temperatures 
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less than 1000 K. For higher temperatures, the two diverge considerably. In the present 
analysis, Monchick's exponential potential is used and, accordingly, equation (C34) has 
been employed for the nitrogen-helium mixture. 


N 2 ~Ar mixture.- Viscosity of nitrogen based on Yos work (ref. 29) is given in the 
preceding section. The viscosity of argon has been obtained from reference 34. A poly- 
nomial fit in temperature is 

M Ar = A + BT + CT 2 + DT 3 + ET 4 + FT 5 — (C37a) 

m 2 

where 


A = 9.25225746 x 10-6 
B = 4.97493659 X 10-8 
C = -7.76309600 X 10" 12 V 
D = 1.70928912 x 10-15 
E = -1.91727863 x 10’ 19 
F = 8.02545298 x 10' 24 ^ 


(C37b) 


The viscosity curve is displayed in figure 21. 

The viscosity of the N 2 -Ar mixture is obtained by using equations (C7) to (C9) along 
with equations (C14) and (C37) The thermal conductivity of the N 2 ~Ar mixture may be 
obtained from equation (C7) with /xj and p 2 replaced by k^ and k 2 obtained from 
equation (C15). The expression for the binary diffusion coefficient (using the Monchick 
exponential potential (ref. 2)) for N 2 ~Ar mixture is 

PD 12 = 1.9816 X10-2 XT 3 / 2 N/sec ( CS B) 

a2l (l,l) 

where T is the temperature in kelvins and 1^ is tabulated as a function of a by 
Monchick. Figure 22 also shows a plot of pD^ as a function of temperature. Poly- 
nomial representation of this function when the collision integral is based on Monchick's 
exponential function is 


pD 12 = A + BT + CT 2 + DT 3 + ET 4 + FT 5 N/sec (300 K i T £ 10 000 K) 
where 

A = -1.42368974 x 10° 

B = 8.50226955 x 10' 3 
C = 1.17660819 x 10“ 5 l 
D = -7.3425527 X 10 -10 
E = 5.01785482 X 10 -14 
F = -1.49447415 X lO' 18 ,, 


(C39a) 


(C39b) 
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Transport by Thermal Diffusion 

Thermal diffusion transport has been included in the analysis of He-N2 mixture. 
The coefficient of thermal diffusion DjT is obtained from the thermal diffusion factor 
through (see ref. 36) 

D^ = pc(l - c)Di2(Q!t) 1 (C40) 

The thermal diffusion factor for the N 2 -He system was computed by using the formula 


(“t)i " ( 6C i2 _ 5 ); 


X 1 S 1 “ X 2 S 2 


(C41) 


+ x 2 2 Q 2 + x 1 x 2 Q 12 

given by Amdur and Mason (ref. 34) based upon Kihara's first approximation. 

The species mole fraction x^ is related to the species mass fraction c^ and the 
species molecular weight by equation (C9). 

The quantities S^, Qj, and Q 12 are given by 

pD 12 


5 «1 2 

S l = 3 = 


4M 1 M 2 A 12 


15M 2 (Mj - M 2 ) 


M i + m 2 11 i H T ( Ml + M 2 ) 2 2 (m j + Mjj" 


(C42) 


Qi = ^r 


M-. 


P D 12 1—2 — 2 8 

A i\/r otv/t ° iv/r tv/t 


( m i 


+ M, 


1 =^ = r — =rr=-( M i + 3M 9 +-M 1 M 9 a: 9 

1 3 /— — \2.ju R'T\ 1 2 5 1 2 12 


(C43) 


Qf2 - 15 l = 


2 

Mj - M 2 \ 32 MjM 2 A 12 . 20 / PD 


vMi + M<- 


+ — Mi Moi 
— — \2 9 1 


12 \ pD 12 


(m x +m 2 )‘ 


p i R t /\ m 2 r T y 


(C44) 


The expressions for S 2 and Q 2 are obtained from those for Sj and by 
interchanging the subscripts, which refer to the molecular species. Conventionally, the 
subscript 1 refers to the heavy molecules (N 2 ) and 2 to the light molecules (He). The 
quantities p.-^ and p 2 (which are approximated by p m ean for N 2 -He mixture) are 
given in equations (C13); pD^ 2 is given in equations (C35). The quantities Cj 2 in equa- 
tion (C41) and A^ 2 in equation (C42) are dependent on the molecular interaction poten- 
tial. These quantities have been evaluated by usin g both MonchickLs_exponential. potential 
(ref. 2) and Lennard -Jones (6-12) potential (ref. 22). In the analysis C^ 2 and A^ 2 , 
based on Monchick's potential, were used for reasons mentioned earlier in this appendix. 
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The polynomial representations of and A^ are when the collision integral 

is based on Monchick's exponential function: 


C* 12 = A + BT + CT 2 + DT 3 + ET 4 + FT 5 

where 

A = 9.4155049 x 10’ 1 
B = -1.35458159 X 10 -5 
C = 3.97227856 x 10" 9 I 
D = -6.93246982 x 10 -13 f 
E = 5.99399177 x 10” 17 
F = -1.99707667 X 10' 21 ,, 

A^ = A + BT + CT 2 + DT 3 + ET 4 + FT 5 

where 

A = 1.1481258 x 10° 

B = 2.63688624 x 10" 5 
C = -7.96613389 x 10" 9 I 
D = 1.39947611 X 10- 12 [ 

E = -1.21357946 X 10‘ 16 
F = 4.04994908 X 10~ 21 > 


(300 K §T 5 10 000 K) (C45a) 


(C45b) 


(300 K s T s 10 000 K) (C45c) 


(C45d) 


and when collision integral is based on Lennard-Jones (6-12) potential: 

C^[ 2 = A + BT + CT 2 + DT 3 + ET 4 + FT 5 (300 K g T g 10 000 K) (C46a) 


where 

A = 9.4356924 x 10' 1 ^ 

B = 6.59447100 X 10" 6 
C = -3.07984512 x 10' 9 I 
D = 6.332799089 x 10' 13 j 
E = -5.91981573 x 10’ 17 
F = 2.05628180 x 10‘ 21 


(C46b) 


Aj 2 = A + BT + CT 2 + DT 3 + ET 4 + FT 5 

where 

A = 1.1024645 x 10° 

B = 3.02127930 x 10" 5 
C = -1.05877112 x 10' 8 I 
D = 1.93950613 x 10' 12 ( 

E = -1.69594890 X 10" 16 
F = 5.62715486 x 10" 21 , 


(300 K g T g 10 000 K) (C46c) 


(C46d) 
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— Lennard-Jones (6-12) potential, ref. 22 

— Monehick' s exponential potential, ref. 2 



Figure 23-- Molecular thermal diffusion factor parameters as functions of temperature. 

These functions are displayed in figure 23. 

The fluid properties evaluated in this appendix may now be used to obtain some 
useful relations required in the solution of the governing equations. 

N 2 ~He mixture.- Employing equations (Cl), (C4), (C7), (C8), (C9), and the assump- 
tion ~ 1^2 ~ M mean > the following relation is obtained for C: 

C = 

p w p w 


p mean 1 1 - 0.8c 1 ~ °’ 857c w 2,282 - 1.282c + 1 - 0.674c 

p w,mean H + 1 1 - 0.8c w 1 - 0.857c c w 1 - c w 

2.282 - 1.282c w + 1 - 0.674c w 


(C47). 


The value of p 


mean 


required in equation (C47) may be obtained from equations (Cl 3). 


Another important quantity required in the analysis is defined by J (the thermal 
-diffusion-factor-): — — — — — — 


D! T i 9T 

T 


(C48) 
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Now, 


1 3T 


1 3H 


0.8 3c 


T 9/3 H + 1 1 - 0.8c 9 /3 


(C49) 


By making use of equations (Cl), (C7), (C8), (C9), (C40), (C49), and the assumption 
Ml ~ M 2 ~ ^mean’ ex P ression (C48) may be written as 


J = 


3.3655 x IQ" 3 c(l - c)(pD 12 )( 


T(7 - 6c) 




1 - c 


mean \2.282 - 1.282c 1 - 0.674c 


x 


_ i 

3H _ ^,2 3c ' 


H + 1 3/3 


C P,2 


Or, by using relation (C4), 


3.3655 X 10-3 


c(l - c)(pD 12 )(Q' T ) 1 

^mean 

(7 - 6c) 

' C | 1 - C \ 

^2.282 - 1.282c 1 - 0.674c/ 


l _1 3H + 0,8 3c \ 

+ 1 3,3 + 1 - 0.8c 3/3 J 


(C50) 


where M mean and pDj 2 ma y be obtained from equations (C13) and (C35). The quantity 
(o'-p)j may be evaluated from equation (C41) by using relations (C9) and (C42) to (C44) 
along with equations (C45). For obtaining T, the following relation, based on defini- 
tions (A21), may be utilized: 

1 - 0.8c w 

T " T » (H + 1} Trnr (C51) 

N 2 -Ar mixture.- Expression similar to equation (C47) may also be obtained for 
the nitrogen -argon mixture by taking into account the appropriate thermodynamic and 
transport properties. The expression is 
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C = -tip 
M W P W 

MjC M 2 (l - c) 

1.188 - 0.188c 1 - 0.168c 

M w,l c w t ^,2^ ~ c w) 

JL .188 - 0.188c w + 1 - 0.168c w 

The quantities and p 2 , required in expression (C52) , may be obtained from equa- 
tions (C14) and (C37). 

Transport due to thermal diffusion was not considered for the N2~Ar mixture and, 
therefore, the thermal diffusion factor J was taken as zero in the analysis for this 


1 1 - 0.5c 


1.429 - 0.429c 


w 


H + 1 1 - 0.5c w 1.429 - 0.429c 


(C52) 
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